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Abstract 

We propose a new method of computing real emission contributions to hard QCD processes. Our 
approach uses sector decomposition of the exclusive final-state phase space to enable extraction of 
all singularities of the real emission matrix elements before integration over any kinematic variable. 
The exact kinematics of the real emission process are preserved in all regions of phase space. 
Traditional approaches to extracting singularities from real emission matrix elements, such as phase 
space slicing and dipole subtraction, require both the determination of counterterms for double 
real emission amplitudes in singular kinematic limits and the integration of these contributions 
analytically to cancel the resulting singularities against virtual corrections. Our method addresses 
both of these issues. The implementation of constraints on the final-state phase space, including 
various jet algorithms, is simple using our approach. We illustrate our method using e + e~ — > jets 
at 0(a 2 s ) as an example. 
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I. INTRODUCTION 



The future high energy collider physics experimental program will measure phenomeno- 
logically interesting quantities with an unprecedented precision. To fully utilize these results, 
accurate theoretical predictions are required. In particular, the large value of the strong cou- 
pling constant as implies that perturbative QCD corrections through next-to-next-to leading 
order (NNLO) in as are needed. The calculation of NNLO QCD corrections has advanced 
rapidly in the past few years. The progress has resulted primarily because of the realization 
that the computation of two-loop virtual corrections can be algorithmically structured and 
automated. These advances culminated in the evaluation of two-loop virtual corrections for 
1 — > 3 and all massless 2 — > 2 scattering processes in perturbative QCD pj. 

Unfortunately, these calculations have not yet produced improved theoretical predictions 
for many observables. The computation of infrared-safe quantities requires two additional 
ingredients: the real-virtual contributions, which denote the one-loop corrections to pro- 
cesses with one additional parton in the final state, and the real-real contributions, which 
denote tree-level processes with two additional partons in the final state. The two com- 
ponents of the real-virtual contributions, one-loop virtual corrections and single emission 
amplitudes, have both been well studied. However, the double emission corrections required 
for the real-real contributions are relatively unknown; the first steps towards understanding 
them have been taken only recently (H, 0] • 

The current state-of-the-art can be illustrated using e + e~ — > 2 jets as an example. Un- 
resolved double real emission corrections appear for the first time at NNLO. The 2 jet 
cross section is currently computed at NNLO by taking the difference of the 0(a s ) e + e~ — ► 
hadrons cross section and the e + e _ — > 3 and 4 jet results at NLO and LO, respectively. Only 
the total 2 jet cross section can be derived using this technique. All information concerning 
the invariant mass and angular distributions of the jets is lost. 

The inability to compute the 2 jet cross section directly at NNLO arises from the poor 
understanding of the singular structure of the double real emission corrections. Infrared and 
collinear singularities cancel between virtual and real corrections only after integration over 
certain kinematic variables makes the 1/e poles in the real emission contributions explicit. 
However, since a primary goal of computing higher order QCD corrections to scattering 
processes is to produce Monte Carlo event generators that correctly describe the kinematics 
of each partonic event, only a restricted region of the final-state phase space can be integrated 
over. Only near the edges of the available phase space, where two or more partons become 
degenerate and combine to form a single jet, can the integration be performed without 
changing the kinematics of the final state. All singularities occur in these limits, and they 
can in principle be extracted and cancelled against those arising from virtual corrections. 
Unfortunately, these singularities overlap; this severly complicates their extraction. 

The existing approaches to computing double real emission corrections extend methods 
used to handle single real emission amplitudes. There are two standard techniques for 
extracting single real emission singularities: phase space slicing jgl [?J and dipole subtrac- 
tion 0, Extending these approaches to double real emission corrections requires two 
non-trivial steps: a determination of the simplified matrix elements that approximate the 
complete double real emission amplitudes in singular kinematic regions, and an integration 
of these matrix elements over the unresolved regions of the multi-particle phase space. The 
difference of the exact and approximate matrix elements is by construction finite, and can be 
integrated numerically. The integration of the approximate matrix elements over the phase 
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space boundaries produces the required 1/e poles that cancel against the virtual corrections. 
Both steps must be completed to obtain an NNLO prediction. Although some progress has 
recently been made 0, E| 0, 11, a functional method for calculating NNLO real emission 
corrections has not yet been demonstrated; a substantial effort is still required to obtain 
phenomenological results. 

We present here a new approach to this problem. We illustrate our technique by consid- 
ering the extraction of singularities from 1 — ^ 4 processes, where the final state particles are 
massless. Our approach is based upon a few ideas. We first derive a factorized parametriza- 
tion of the four particle phase space following a simple procedure. We, then, use sector 
decomposition [ll|, [l^, *° separate the overlapping divergences which appear in the 
double real emission matrix elements. After this separation is performed, the phase space 
singularities can be extracted using standard expansions in terms of plus distributions. The 
processes of finding the required sectors and extracting the singularities are completely au- 
tomated. The resulting matrix elements are finite and fully differential, and can be used 
to create NNLO Monte Carlo event generators. We discuss in some detail the example of 
e + e~ — > hadrons at C?(a|), which includes the 2 jet cross section at NNLO, the three jet 
cross section at NLO, and the 4 jet cross section at LO. 

The paper is organized as follows. In Section II we introduce our method by considering 
e + e~ — > 2 jets at NLO and e + e~ — > 3 jets at tree level. We begin our discussion of the 2 jet 
cross section at NNLO in Section III by describing our parameterization of the four-particle 
phase space. We also explain how we use sector decomposition to separate the overlapping 
singularities that appear in the matrix elements. In Section IV we apply our technique to the 
two most difficult interferences that appear in the double real emission contributions. After 
demonstrating that our method is powerful enough to handle the most complicated scenario, 
we apply it to a simple but realistic example in Section V: the Nf dependent contribution 
to e + e~ — > hadrons at 0(a 2 s ). This process contributes to the 2 jet cross section at NNLO, 
the 3 jet cross section at NLO, and the 4 jet cross section at LO. Finally, we present our 
conclusions and discuss future prospects in Section VI. 



II. THE NLO EXAMPLE 

We begin by considering the 0(a s ) contribution to e + e~ — > hadrons, which contains both 
the NLO correction to e + e _ — > 2 jets and the LO contribution to e + e~ — > 3 jets. Although 
many of the complexities of the NNLO case are absent in this calculation, it illustrates 
several important features of our method. At the partonic level, we must compute the 
one-loop virtual corrections to e + e~ — > qq and the real emission process e + e~ — > qqg. 

We consider first the real emission correction e + e~ — > qqg. The kinematics of the final 
state is fully described by the invariant masses s qq -, s qg , and Sg g , which satisfy the constraint 

s qq + s qg + s qg = S. (1) 

Here, s is the center of mass energy squared of the colliding electron and positron. The 
three particle phase-space can be written as 

[dq] [dq] [dg] (2n) d 5^ ( Pl + p 2 - q - q - g) = 

1 ^T^T) / dA i d W(l " Ai)- e A 2 -(l - A 2 ) 1_2e , (2) 



o 
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where d = 4 — 2e, [dk] = d d 1 k/{2ir) d 1 fc , TI2 is the integrated phase-space of the two 
massless particles, 

><-2 



and Qd is the solid angle in d dimensions, 



2ni 



(3) 



(4) 



The invariant masses have the following expressions in terms of Ai and A2: 



>qq 



s(l-A 2 )(l-Ai), 



S(l — A2)Al, Sgg — SA 2 - 



(5) 



>qg °qg 

In what follows we set s — 1 for simplicity, and restore the correct dimensions in final results. 

The matrix element for the 7* — > qqg process is given by two diagrams. Upon squaring 
these and using the expressions for the invariant masses given in Eq. El we derive 
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2(1 - Ai)(l - A 2 ) + Xj + A?(l - A 2 ) 2 - e [Ax + A 2 - A X A 2 ] 



(6) 



AiA 2 (l — A 2 

After substituting the expression for the matrix element squared into the three-particle 
phase-space, we arrive at the expression 



Tin 



'ha ha dV * _J 

12 dA x dA 2 (4vr) rf / 2 T(l - e) 



dAidA 2 (AiA 2 ) 



-e-l, 



l-A 1 )- e (l-A 2 )-^(A 1 ,A 



(7) 

where g>(Ai,A 2 ) is a non-singular function of the A». The phase space singularities in the 
above expression can be extracted before integration by using the standard decomposition 
in terms of plus distributions: 
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where a plus distribution is defined via 



dX 



ln w (A) 
A 



/(A) 



dAln n (A) 
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Substituting this decomposition into Eq. [7J we derive the following expression for the real 
emission cross section: 

d 2 a R 64vra s T(l + e) f 5(X 1 )8(X 2 ) , 1 [ 5(Xi 



dAidA 2 
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5(A 2 ) + |A 2 - 1 
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(2 - 2A 2 + Aj) ln(l - A 2 
A 2 
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8{X 1 )8{X 2 
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o"o is the tree level cross section for e + e~ — > qq: cr = AnaEMQl/ s. For the calculation of 
the NLO corrections to the 2jet cross section we also require the virtual corrections to the 
e + e" — > qq process; we find |6[ 

da v 64 TO r_(l + e)| 1 1 2tt 2 5\ 



dA 1 dA 2 3 u (4tt) o! / 2 [ e 2 2e 3 2 

We now discuss the calculation of the n jet cross section; here, n equals either 2 or 3. We 
introduce the jet function 

F^\s qq , s qg , s qg ) = Fj n) ((l - Ai)(l - A 2 ), Ax(l - A 2 ), A 2 ). (12) 

The n jet cross section becomes 



(n) 



X^fPm L M x 2) + + . (13 ) 



It is clear from the expressions in Eqs. an d ^2 that the 1/e 2 poles cancel when ay and 
cxr are combined. The 1/e poles in ctr require that either Ai or A 2 vanish. The jet function 
becomes either Fj (s qq -, s qg ,0) or Fj(s qq , 0, s qg ) in these cases, i.e., a 2 jet configuration 
is always obtained. The 1/e poles of both <jy and ctr occur in the 2 jet cross section; they 
cancel after integrating over Ai and A 2 , as required for infrared safe observables. Dropping 
the poles in e, we can write the n jet cross section as an integral over the finite component 
of the partonic cross sections: 



/ dA x dA 2 

"^J @ finite- 

(14) 

Jo 



It is straightforward to check that this correctly reproduces known 2 and 3 jet cross sections 
for standard jet functions 

Several important aspects of this result generalize immediately to NNLO calculations. 
We were able to extract the singularities in e without performing any integrations. The 
cancellation of the poles in e can be checked numerically, and these terms can then be 
discarded. We note that by casting the subtraction operations needed for extracting the 
e poles in terms of plus distributions, we have gained the flexibility to combine our result 
with any jet function and with any constraint on an infrared safe differential quantity. This 
greatly simplifies the calculation of NNLO cross sections. 



III. FOUR PARTICLE PHASE-SPACE: PARAMETERIZATION AND SECTOR 
DECOMPOSITION 

We present here our parameterization of the four particle phase space. After deriving 
the relevant formulae, we discuss the complications that arise when we attempt to extract 
phase space singularities using the method discussed in the previous Section. We show how 
sector decomposition of the phase space solves these problems, and apply the technique to 
a few examples. 

We begin with the following expression for the four particle phase space: 

J 4 = / [dp 1 }[dp2}[dp 3 }[dp 4 }(27r) d 5^(p-p 1 - p 2 - p 3 -p 4 ). (15) 



5 



This phase space is described by five independent invariant masses. A convenient set is 
{S134, s 2 34, S34, si 3 , S24}, where Sy = (p, + p.,) 2 and s ijk = (pi + pj + p k f . We split the above 
integral into three sub-integrals: 



-1) 



3d-4 



dS234dS34dSi34dS23dSi 3 Ji/2^3, 



where 



and 



h = j d d Pl d d Q 2M 5^ (p- Pl - Q234) S (p 2 ) 5 (Q 2 34 - s 2U ) , 
' d d p 2 d d Q M 5^ (g 234 -P2- Qu) 5 (pi) 5 (Q 2 4 - s u ) , 

h = J d d p 3 d d P4 5^ (Q 34 -ps- p 4 ) 5 (pf) 5 (pi) . 



(16) 

(17) 
(18) 

(19) 



We now constrain the integrations over the by introducing the following delta functions 
into Z 4 : 

S (s 2 34 - Q234) $ ( s 34 - Qti) & (si34 - Qk) 5 ( S 23 ~ 2p 2 • Pa) & («13 - 2f>i • p 3 ) • (20) 

To derive representations of these integrals from which the phase space singularities can 
be conveniently extracted, we bring the limits of integration for each integral from to 1 
using transformations of the form Sy = \ (sfj — + s~-, where denote the maximum 
and minimum values of the corresponding invariant masses. Using the delta functions to 
simplify the integrations, performing the transformation to the variables Aj, and including 
the Jacobian \dsij/d\k\, we arrive at 

J 4 = Na C dA 1 dA 2 dA 3 dA4dA5 5(Ai - \[)5{\ 2 - A' 2 )5(A 3 - A 3 )<5(A 4 - Al)5(A 5 - A' 5 ) 
Jo 

x [A x (1 - Ai) (1 - A 2 )] x - 2e [A 2 A 3 (1 - A3) A 4 (1 - A 4 )P [A 5 (1 - A 5 )]- 1/2 - e ■ (21) 
We have extracted the overall normalization 
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(2tt) M - 2 r 2 (l + e)r 4 (l-e)r 2 (l/2-e) 



(22) 



The invariant masses have the following expressions in terms of the A, (with s = 1): 

•S234 = Ai, 

S34 = A1A2, 

S23 — Ai (1 — A2) A 4 , 

S134 = A 2 + A 3 (1 - Ai) (1 - A 2 ) , 



Sl3 



A 5 (i 



13 



'13 



+ S 



13' 



with 



'13 



Ai) 



A 3 A 4 + A 2 (1 - A 3 ) (1 - A 4 ) ± 2v%A 3 (1 - A 3 ) A 4 (1 - A 4 ) 



(23) 



(24) 
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Difficulties arise when we substitute the matrix elements into the four particle phase 
space and attempt to expand the expression using Eq. |H1 We will discuss in Section V 
the Nf contributions to e + e~ — > 2 jets; the matrix elements for the double real emission 
contribution to this process contain denominators of the form 1/53432343134. Using Eq. |22l 
this becomes 

1 = 1 . (25) 

S34S234S134 A?A 2 [A 2 + A 3 (1-A 1 )(1-A 2 )]' v ; 

The third term in this denominator is singular when e.g. both A 2 , A3 — > 0, but not when 
only one does. If we combine the denominator with the integration measure in Eq. |2~T1 
and attempt to naively expand X^ 1 ^ 6 — > — <5(A 2 )/e + AjT 6 — > 1 — eln(A 3 ) + we will 
find unregulated singularities as A3 — > 0. The most convenient method for separating the 
overlapping singularities in A 2 and A3 is sector decomposition 11, I2,[l3j]. To illustrate this 
technique, a simple example suffices. We consider the integral 

1 

dx dy x~ 1 ~ e y~ 1 ~ e (x + y) e . (26) 



The 1/x and 1/y factors cannot be expanded in plus distributions, as the logarithms from 
the expansion of x + y will produce singular terms. We split this integral into two parts, 

h = f 1 dx f dy x~ l - e y- 1 ^ (x + y)-\ h= C dy f dx x^y' 1 ^ (x + y)~ e . (27) 
Jo Jo Jo Jo 

111 Ii we set y = y/x, and in I 2 we set x = x/y. Performing these variable changes, we find 

h = f 1 dx dy x- 1 - 3 ^- 1 -" (I + y y e , I 2 = f 1 dxdyy- 1 ' 3 ^- 1 - 6 (1 + x)~ e . (28) 
Jo Jo 

The singularities in x and y are now separated in each integral, and can be extracted using 
Eq. El 

One great advantage of this technique is the ease with which it can be automated. The 
rules to determine when a term requires sector decomposition are simple; if the expression 
becomes singular when two (or more) variables x, y — > 0, but remains finite when either 
x —>■ or y — > 0, then the transformations discussed below Eq. E?| should be performed. 
Another advantage of sector decomposition is that it can be applied to fractional powers in 
addition to denominators, as illustrated in the example above. 

We now discuss the application of this method to the denominator in Eq. to show 
how it works in practice. It is convenient to first separate the two singularities that can 
occur if x — > or x —>■ 1 by splitting the integration as 

r l f l/2 f l 

dx^ dx+ dx , (29) 

Jo Jo Jl/2 

and changing x — > x in the second integration so that x — 1 is mapped to x =0. Doing so 
for the three variables Ai, A 2 , and A 3 produces eight sectors. We focus on the sector where 
originally Ai, A 2 , A3 < 1/2. The denominator has the form 



V 



A?A 2 A 2 + A 3 (l-^)(l-4f) 



(30) 
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We now perform a sector decomposition in the variables A2 and A3, using the transformations 
given below equation Eq. in sector a we set A2 — > A2A3, and in sector b we set A3 



A3 A 



3^2- 



Combining the Jacobian of the variable change with the denominator, we find the following 
expressions in each sector: 



V a 



1 



Af A2A3 



A 2 + 



A2A3 - 

2 - 



1 



A?A 2 1 + A 3 (1-^)(1- 



A2 

2 • 



(31) 



The terms in brackets are now finite in all limits; the denominators can be combined with 
the phase space measure, and the standard decomposition in plus distributions can be used 
to extract singularities. Note that the above transformations must also be performed in the 
integration measure. After splitting the integration as in Eq. the measure contains terms 
of the form (1 — Aj/2). After sector decomposition, these become (1 — AjAy/2). If we had 
not performed this split, we would have produced terms of the form (1 — AjAj). These are 
potentially singular when Aj, Xj —>■ 1, and would require further sector decomposition. 

We must discuss two subtleties that can occur when using the method presented above. 
The representation of Z 4 we have derived is convenient for expressions that do not contain 
S13 or s 14 (see Eq.l23|) in the denominator. In such terms, it is difficult to extract singularities 
in A 2 , A 3 , A 4 that appear after integrating over A 5 . One can always re- map the momenta 
of the final state particles in a given diagram in such a way that si 4 never appears in the 
denominator. For those terms that contain S13, we first bring the limits of the S13 integration 
from to 1 using the transformation 



A, 



$13 — $13 $13 



s 13 ~ s 13 S 13 

We then derive the following expression for the four particle phase space: 
-I 
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(32) 



J 4 = A/" 4 f dA 1 dA 2 dA 3 dA 4 dA 5 5(Ai - A' X )5(A 2 - A 2 )5(A 3 - A 3 )5(A 4 - A 4 )5(A 5 - A 5 



x [Ai (1 - Ax) (1 - \ 2 )} 1 - 2e [A2A3 (1 - A 3 ) A 4 (1 - A 4 )]" 



A, 1 - A, 



-1/2-e 



x si 3 (A E 



S 13 S 13 



-1/2-e 



13 



'13 



'13/ 



2c 



(33) 



Factors of S13 that appear in the denominator are cancelled by the Jacobian of the non-linear 
transformation of Eq. |H2l and the remaining A5 integration does not produce dangerous 
singularities. Therefore, sector decomposition of the A5 integration is never required. 
One further complication exists. Using the expressions in Eq. we find 



S 13 S 13 



-1/2-e 



(1 - Ax)- 1 ^ |A 2 (1 - A 3 ) (1 - A 4 ) - A 3 A 4 



-l-2e 



(34) 



This expression is singular on a manifold of points in the interior of the phase space. We 
wish to move these singularities to the boundary of the integration region. To do so, we first 
note that the singularity occurs when 



A, 



A^ 



A 2 (l-A 3 ) 
A3 + A2 (1 — A3 



(35) 



S 



this value is always in the integration region. We can therefore split the A4 integration into 
two parts, 

/ dA 4 = / dA 4 + / dA 4 , (36) 
Jo Jo Jx\ 

and then bring the integration limits back to and 1. Doing so produces two integrals, T\ 
and Z4, with all singularities moved to the boundaries of the integration regions; these can 
be extracted using the sector decomposition technique discussed above. 

To summarize, we derive analytic results for the double-real radiation corrections by 
following these steps: 

• We derive a factorized parametrization of the 1 — > 4 phase-space as in Eq.l2Tlor Eq.lBlfl 
in terms of kinematic variables which range from to 1 . 

• We remove singularities from inside the allowed phase-space region to the boundaries 
by splitting appropriately the integrations and mapping them back to the [0, 1] interval. 

• We then apply sector decomposition to disentangle the overlapping singularities. 

• Finally, we extract the e poles in terms of plus distributions. 



IV. THE UNITARITY CHECK 

Having discussed the four particle phase space and the technique of sector decomposition, 
we now illustrate our method by considering two examples of double real emission integrals 
with four propagators. These are the most complicated phase space integrals that appear in 
1 — > 4 processes. We check our calculation of the double real emission corrections using their 
contributions to the imaginary parts of three-loop propagator diagrams. From the optical 
theorem we know that the imaginary parts of such diagrams are given by the sum of all cuts, 
where all possible combinations of internal propagators are put on-shell. The required cuts 
also include real-virtual and virtual-virtual ones. These are simple to compute, as are the 
imaginary parts of the propagator diagrams. We can therefore derive analytic expressions 
for the real-real cuts, which we can compare with the results we obtain using the methods 
presented in the previous Sections. The checks we perform in this section involve inclusive 
integrations over the real emission phase space, in order to compare with the imaginary 
part of the relevant propagator diagrams. We will demonstrate in the next Section that the 
implementation of jet functions into our method is simple in both principle and practice. 

We begin by considering the maximally planar propagator integrals which arise from 
Feynman graphs as the one shown in Fig. 1. The sum of all cuts for these diagrams con- 
tributes to the e + e~ total cross-section. Here, we examine the underlying scalar integral of 




FIG. 1: A planar propagator-type diagram contributing to the e + e cross-section at NNLO. 
these diagrams after we set their numerator to one. The analytic expression for this integral 
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can be found using MINCER [151 ]: it is finite, and therefore its imaginary part vanishes. The 
sum of all possible cuts of this diagram must also vanish. This diagram has three distinct 
virtual-virtual cuts: the interference of a two-loop planar vertex correction to 7* — > qq with 
the tree level contribution diagram together with its complex conjugate, and the square of 
the one-loop vertex correction to this process. It has four real-virtual cuts: two copies of 
the tree-level contribution to 7* — > qqg interfered with the vertex correction to this process 
with the gluon radiated off the opposite quark leg, and two copies of its complex conjugate. 
Finally, it has two real-real cuts: the tree-level contribution to 7* — > qqgg with both gluons 
radiated from a single quark line interfered with the diagram where both gluons are radiated 
from the opposite quark line, together with the complex conjugate of this contribution. 
The first virtual-virtual term, the two-loop vertex correction and its complex conjugate, 

is ET 



V - ( T{1 + e) Y<!l ( 1 4C2 + 10Cs 38C2 ^ (37) 

where 1Z 2 is the the two-particle phase space introduced in Section II. The second virtual- 
virtual contribution, the square of the one-loop vertex correction, is 

The four real-virtual contributions give 

The sum of the above contributions with a minus sign should equal the sum of the real-real 
cuts of this diagram. We obtain 



R 



/£(l + e)y fj_ _ 2<2 
^ (4vr)<V 2 J 2 \2e* e 2 

/r(l + e)\ 2 _ /0.5 3.2899 




n A — ~ + — + 107 - 15 - ) • ( 4 °) 



V (47r) d / 2 J z V e 4 e 

In terms of the invariant masses introduced in the previous Section, this contribution should 
be equal to 

iC m = 2 ( ), (41) 

where the factor of two indicates the sum of both real-real cuts. This integral involves the in- 
variant mass S\3] its calculation therefore requires extensive use of sector decomposition and 
the reparameterization A5 — > A5 discussed in the previous Section. Using these techniques, 
and numerically integrating the result using VEGAS 17], we obtain 

_/r(l + e)\ 2 /0.5 (-0.8 ±9.2) • 1(T 5 3.2909 ± 0.0018 16.827 ±0.010 

R r> -{-J^W) n * {-^ + i5 = + 



6 2 



+107.12 ±0.07). (42) 
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We have included the VEGAS errors for those terms which require non-trivial integrations. 
The agreement between the analytic and the numerical results is better than 0.1% for all 
terms considered, and the differences are consistent with the integration errors. 

We now consider the non-planar topologies shown in Fig. 2. Once again we focus on the 
scalar integral which is obtained from these topologies by setting their numerator to unity. 
The imaginary part of this integral vanishes; therefore, we can again verify our calculation of 
the double real emission contribution using sector decomposition by checking the cancellation 
of all possible cuts. This diagram has two virtual-virtual cuts: the interference of the two- 
loop non-planar vertex correction with the tree level 7* — > qq diagram, and its complex 
conjugate. It has four real-virtual cuts, each of which involves the one-loop box correction 
to 7* — > qqg interfered with the tree-level contribution. Finally, it has five real-real cuts, of 
three distinct types: two cuts involving the emission of two gluons off a single quark line 
interfered with the emission of two gluons off the opposite quark line; two cuts where a 
radiated gluon splits into a qq pair; one cut where a gluon is radiated from each quark line. 
Since we are considering only scalar diagrams, the first four cuts give identical answers. 




FIG. 2: A non-planar propagator- type diagram contributing to the e + e cross-section at NNLO. 
The virtual-virtual cuts sum to E3 



V, 



np 



£(1 + 6) 

( 47r )d/2 



38C 2 



54(3 + 966CJ 



5 



(43) 



The real-virtual contribution involves the integration of the one-loop box diagram with one 
leg off-shell. There are two distinct ways of computing this diagram: either analytically 
using e.g. a Mellin-Barnes representation, or numerically using sector decomposition in a 
fashion identical to our approach to the real-real terms. The second method is particularly 
convenient, as it allows restrictions on the final-state phase space to be imposed easily. We 
now illustrate this technique. 

The expression for the one-loop scalar box diagram with one external leg off-shell, valid 
to all orders in e, is 



B 



r(l-e) 2 r(l + e) 1 



e 2( 47r )d/2 r(l-2e) 
+ (- S )- e F 21 (l,-e,l-e, 



st 



e, 



■-) - (-M 2 
t 



-t)" e F 21 (l,-e,l 

F 21 (l,-e,l-e, 



'~sT' 



(44) 



Here, M 2 is the virtuality of the off-shell leg, and s, t, u are the usual Mandelstam variables. 



We set M 2 = 1 in what follows and choose the Mandelstam variables to be s = s 



A, 



t = Sq g = (1 — A 2 )Ai, u = Sqg = (1 — A 2 )(l — Ai), where we have used our notation from Section 



II. The expression in Eq. 1441 must be integrated over the three particle phase, together with 
an additional propagator l/sq g coming from the interference with the tree- level diagram, 
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and the three-particle phase space in Eq. EJ It is convenient to proceed in the following way. 
Using the integral representation for the hypergeometric function, 



T(i-e) } , r- 1 



Ml,-e,l-e,z) = ^yd*— , (45) 

v ' o 

we can write the required integrals over the phase space and the auxiliary variable t in a form 
that is directly amenable to sector decomposition. It is then a simple task to numerically 
compute the expansion of the real-virtual corrections in powers of e using the techniques 
described above. The analytic result for the real-virtual cut was derived for the purpose of 
checking our calculation based on sector decomposition; summing the four cuts, we obtain 

TtV ( V{lJre) W ( 10 128C 2 356C 3 1108£\ 
RVnp -{-{^W) n2 {~^ + — + — 5~ J' (46) 

Requiring the cancellation of all cuts, we derive the following analytic expression for the 
real-real contribution: 

2 /o nnA onoA , AOf2\ 



£(l + e)r 8 90C 2 302C 3 142& 

Rnp -{-J^WJ n2 \7~—i ~ + — 

2 



/T(l + e)\ * (8 148.04 363.02 n . 

K 2 1 5 ± 76.84 . (47) 



" ^ (4tt)<V 2 ; Ve 4 e 2 
This result should be equal to the following sum of real-real integrals: 

RZ m = ( ) +4/ ) . (48) 

Computing the above integrals using sector decomposition, we obtain 

/T(l + e)\ 2 ^ [8 (5.8 ±6.4) • 10~ 4 148.06 ± 0.03 363.03 ±0.07 



R 



\ (4vr) d / 2 J z \e 4 e 3 



e 



2 



±77.1 ± 0.4 ) . (49) 



The result found using sector decomposition is again consistent with that found by demand- 
ing the cancellation of all cuts, although the numerical precision of the finite piece is slightly 
worse. This can be improved with a more sophisticated numerical integration technique. 

We therefore conclude that our method can accomodate the most difficult real-real phase 
space integrals needed for 1 — > 4 processes. We next consider in detail the Nf dependent 
contributions to the e + e~ — > 2, 3, and 4 jet cross sections. This example addresses the two 
remaining issues we must confront to fully validate our method: that jet functions can be 
implemented simply, and that bookkeeping of the sectors (i.e., the expressions for the s»j in 
each sector as a function of the rescaled A») can be peformed. 



V. N f DEPENDENT CONTRIBUTION TO e + e~ -> 2,3,4 jets 

In this Section we illustrate our method in a realistic NNLO example. We compute the 
Nf dependent contributions to the 2 jet cross section at NNLO. When we wish to compute 
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jet cross sections, we must include in our matrix elements a jet function, denoted by Fj 
in Section II, that determines whether a given configuration contains 2, 3, or 4 jets. This 
function takes the invariant masses of all partonic pairs as its arguments. After splitting 
our result for the double real emission contribution into sectors, the invariant masses take 
different forms in terms of the A; in each sector. This presents bookkeeping issues that must 
be addressed. We prove here that we can handle this problem by considering a realistic 
example. Taken in conjunction with the calculation in the previous Section of the most 
difficult integrals that appear in e + e~ — ► 2 jets, the reader should be convinced of the power 
of our approach. 

An example of the diagrams that contribute to the Nf dependent terms of e + e~ — ► jets 
at 0(a 2 s ) is shown in Fig. 3. This diagram, together with the remaining contributions where 
the internal bubble is attached to a single quark line, contains both virtual- virtual and real- 
real cuts. At O(o.g) we must also consider the coupling constant renormalization of the 
0(as) result. To present numerical results we must also choose a jet algorithm; we use the 
JADE algorithm, with a separation parameter y = 0.1. 




FIG. 3: An example of an Nf dependent diagram at NNLO. 

The virtual correction only contributes to the two-jet configuration. Therefore, we write 

°v c » T (<*>\ 2 ( 1 11 l ( 11 > 269\ 5423 13 , 121 > \ 

V a = S '* N ' IT) (l8? + 54? + 1 ("I8 C2 + 32i) + 1944 " 27 6 " M*) ' (50) 

where the Kronecker delta indicates the restriction to the 2 jet cross section. We compute 
the double real emission using the approach described in the previous sections, and extract 
the singularities prior to integration over any kinematic variables. Implementing the jet 
algorithm and performing the integrations over the five-dimensional phase-space numerically, 
we obtain 



aR Nf (—S* (8 



(5.5553 ± 0.0005) • 10~ 2 0.20369 ± 0.00005 0.4180 ± 0.0005 



+4.808 ± 0.003 ] - S ji3 (0-41005^0.00016 + ^ ± ^ 

+ (1.4561 ± 0.0018) • 10^-4 ) . (51) 

The 0(a s ) cross-section combines the virtual correction and the single real emission. We 
need this contribution to 0(e) to derive its contribution to the NNLO cross section. Using 
the results in Section II, we derive 

/ — W/ ((-1.4597 ±0.0013- (9.242 ±0.004) e) 8 j>2 + 



a \rr 



(2.4575 ±0.0012 + (6.115 ±0.003) e) 5 ji3 ) . (52) 
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The 0(a 2 ) contribution to the e + e annihilation into hadrons is then written as 



a 



(2) 



OV + or - 



A) fat 



71 



a 



(i) 



(53) 



The last term comes from the renormalization of the strong coupling constant in the 0(a s ) 
cross section; we need only keep the —Nf/6 term in the beta-function. Adding these con- 
tributions, we obtain 



a 



(2) 



00 



N f 



(I; 



J.2 



[2.Q ±4.6) • 1(T 6 (1.4 ±5.5) • 10~ 5 (3.1 ± 5.2) • 10" 



±1.799 ± 0.003 ) ± 6 jt3 i 

± (1.456 ± 0.002) • lO~ 3 5 jA } . 



-0.5 ±2.6) ■ 10' 



1.917± 0.017 



(54) 



As we see, the divergences associated with various pieces disappear, with small remnants 
consistent with the integration errors. The cancellation occurs independently for the 2 and 
3 jet cross-sections, as required. Finally, adding together the 2, 3, and 4 jet cross sections, 
we obtain the total hadronic cross section 



a 



(2) 



00 



-0.117 ±0.003) N f 



a. 



71 



(55) 



which agrees with the known analytic result 

r( 2 ) /2 _ 11 



(J 



3C3 



-0.115 Nf ( — 

7T 



(56) 



Again, the integration error of the finite piece can be improved with a more sophisticated 
numerical integration technique. We conclude that our method can be applied succesfully 
to compute differential quantities. 



VI. CONCLUSIONS 

We have presented a new technique for computing double real emission corrections at 
NNLO. Our method uses sector decomposition of the four particle phase space, together 
with an expansion in plus distributions, to extract the phase-space singularities without 
any analytic integrations, and preserves the exact kinematics of the partonic event. The 
expressions for the matrix elements obtained with this approach can be used as building 
blocks for Monte Carlo event generators. 

A phenomenologically attractive feature of our method is that constraints on the final- 
state phase space, including various jet algorithms, can be implemented simply. This makes 
it possible to study radiative corrections to quantities of direct experimental relevance. The 
method is completely automated, and flexible. It can be applied to any QCD or electroweak 
process with massless particles in the final state, where the singularities from double-real 
unresolved radiation must be extracted. 

We have illustrated our approach using e + e~ — ► jets at 0{a 2 s ) as an example. We have 
considered the most complicated phase space integrals that appear. We have explicitly 
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checked our results for those integrals by performing an inclusive numerical integration 
over the phase-space and comparing with analyting results obtained by using the unitarity 
method. We have also demonstrated that our method is capable of calculating differential 
quantities at NNLO by deriving the Nf dependent contributions to the 0(a 2 s ) cross section 
for e + e~ — > jets; this includes e + e~ — > 2 jets at NNLO, e + e~ — > 3 jets at NLO, and e + e _ — > 
4 jets at LO. The 1/e poles were cancelled numerically, and the finite piece for the inclusive 
cross section agrees with results in the literature. 

Results for the non-Ay contributions will be given elsewhere. We have already presented 
here the calculation of the most difficult contributions needed for these terms. In addition, 
the bookkeeping of the various sector decompositions, and the numerical integrations, have 
already been addressed here. 

While a direct application of our formalism to more complicated phase-spaces is a viable 
option, one could also use it profitably in conjunction with a dipole formalism. Sector 
decomposition can be applied to the dipoles to extract the 1/e singularities they contain, 
without the need for an analytic integration. The remainder can be integrated numerically. 
Sector decomposition of the finite terms should also improve the numerical stability of the 
dipole approach. 

There are several possibilities to develop the method further. It is interesting to investi- 
gate its direct application to 1 — > 5 processes. It is also important for many applications to 
study the factorization of the phase space when massive particles appear in the final state. 
Although the parameterization of the phase space is certainly more complicated in those 
cases, we do not anticipate any significant limitations of the method. We also expect that 
the number of required sector decompositions will be reduced in the presence of massive 
particles. 

Our method is a promising new technique for computing real radiation contributions to 
NNLO cross sections, and allows to obtain phenomenological results vital for the future of 
precision high energy physics. We look forward to the application of our method to many 
important collider physics processes. 
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Note Added: While this paper was being completed, a new paper appeared which discusses 
the application of sector decomposition to inclusive phase space integrals (l^ |. 
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